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Abstract. Sum rule constraints dominated by the independent high-scale input, a s (Mz), are shown 
to be satisfied by I = 1 spectral data from hadronic T decays, but violated by the pre-2005 electropro¬ 
duction (EM) cross-section data. Determinations of the Standard Model (SM) hadronic contribution 
to (g — 2 )jj incorporating % decay data are thus favored over those based solely on EM data, implying 
a SM prediction for (g — 2) u in agreement with current experimental results. 

After accounting for the accurately known, purely leptonic contributions, the largest 
of the SM contributions to = (g — 2)^/2 is that due to leading (LO) hadronic vacuum 

polarization (VP), This contribution can be expressed as a weighted integral, 

with known kernel, over the electromagnetic (EM) spectral function Pem{s). The un¬ 
certainty on \< 2 n\ had which results is at the ~ 1% level (comparable to the 0.5 ppm 
experimental uncertainty on a^ [ 1]) and dominates the uncertainty in the SM prediction 
for afi Hi]. This determination can, in principle, be improved using CVC, together with 
data on the I — 1 spectral function, p/=i(s), obtained from hadronic T decay. (Isospin¬ 
breaking (IB) corrections Hi are needed for accuracy at the 1% level.) Unfortunately, 
the IB-corrected r data is in significant disagreement with the corresponding 7=1 EM 
data, p/'ii'j lying uniformly below pj =i in both the 4k region and that part of the 2 n region 

above the p but below 1 GeV 2 fl. The determination of employing only EM 

data leads to a SM a^ prediction ~ 2.5<7 below experiment, while the alternate determi¬ 
nation incorporating x decay data is compatible with experiment at the ~ 1 o level 0] • 
We show here that sum rule constraints strongly favor the T-based determination. 

For vv(.v) any function analytic in the region ,v < M, with M > .sp, the 7 = 1 vector 
(V) current and EM correlators, 11(7), and corresponding spectral functions, p ( 7 ), satisfy 
the FESR relations fg°w(s)p(s) ds = -if|,|= io vv ( 5 ) n ( 7 ) ds. To suppress duality vio¬ 
lation and allow the use of the OPE on the RHS, w(7) should satisfy w(s — 7o) = 0 0]. At 
scales of ~ a few GeV 2 , the OPE representation for V current correlators is essentially 
entirely dominated by its leading D = 0 component, and hence determined by the single 
input parameter a s {Mz), whose value is known from independent high-scale studies [7]. 
If one works with weights w(y), y — s/sq, which satisfy w(l) = 0, and are non-negative 
and monotonically decreasing for 0 < y < 1, the fact that the EM version of p/=i lies 
uniformly below the corresponding x version implies that 

• the normalization and slope with respect to 7o of the EM-based spectral integrals, 
for a given w(y), should be too low relative to OPE predictions if it is the T-based 


TABLE 1 . a s {Mz ) from fits to the .vq = m~ experimental EM and T 
spectral integrals, with central D = 2,4 OPE input values 


Weight 

[ a s (Mz)\ EM 


W 1 

n 11 qq +0.0030 

O.11a8_ 000 3 5 

0 1212 +00027 
u.iz,iz_ 00032 

W6 

0 1150+ 0 0022 
U.11JU_ 00026 

0 1 1 +0-0020 
U.11^J_ 00022 


data which is correct and 

• the normalization and slope with respect to so of the T-based spectral integrals for a 
given w(y) should be similarly too high if it is the EM-based data which is correct. 

Sum rule tests of the two data sets have been performed for a number of different 
w(y). See Ref. H] for details on the OPE and spectral integral inputs (including short- 
and long-distance IB corrections for the r data and input relevant to the small D > 0 
OPE contributions). It is found that both the normalization and slope with respect to sq 
of the T based spectral integrals are in excellent agreement with OPE expectations, while 
both the normalization and slope of the EM-based spectral integrals are low, particularly 
if one uses the pre-2005 EM nn spectral data. 

Table [T]quantifies the EM normalization problem, giving the values of a s (Mz) needed 
to bring the OPE and spectral integrals into agreement for the T and EM cases, at the 
maximum common available scale, sq — m Results are shown for w(y) — vtq (y) = 1 — y 
and W(,{y) = 1 — 6y/5 +y 6 /5 (which have zeros of order 1, 2, respectively, at s = so). 
The entries are to be compared, e.g., to the high-scale average, a s (Mz) = 0.1195 ± 
0.0016, obtained by excluding (i) z decay input and (ii) the erroneous heavy quarkonium 
input Dafl from the PDG04 average. Table [2] similarly quantifies the EM slope problem. 
In the OPE column, “indep” and “fit” label results obtained using (i) the independent 
high-scale a s (Mz) input and (ii) the fitted a s (Mz) values from Table Q] respectively. 
We see that lowering the input a s (M z ) value to accommodate the EM spectral integral 
normalizations at ,vo ~ m\ does not resolve the slope problem for the EM spectral 
integrals. 

For illustration purposes, the results for the w(, EM case are also displayed in Fig. [H 
The dotted and solid lines give the central OPE and OPE error bounds, respectively, the 
solid dots (with error bars) the EM spectral integrals, with pre-2005 nn input. The open 
circles show the shifted “EM” spectral integrals obtained by replacing the EM In and 
An data with the equivalent z data. We see that both the EM slope and normalization 
problems are resolved if, where the I — IV part of the EM and t data disagree, the z 
data are taken to be correct. The recent 2005 SND EM nn results are compatible 
with the corresponding z results, further strengthening the case for the reliability of the 
T data, and the agreement between the SM prediction and experimental result for a a ■ 
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TABLE 2. Slopes wrt so of the EM OPE and spectral integrals 


Weight 

S exp 

cc s (M z ) 

SoPE 

W 1 

.00872 ± .00026 

indep 

.00943 ± .00008 



fit 

.00934 ±.00008 

W6 

.00762 ±.00017 

indep 

.00811 ±.00009 



fit 

.00805 ± .00009 
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FIGURE 1. OPE vs. spectral integrals for the wg FESR. All notation as described in the text. 
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